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ABSTRACT 

. We demonstrate that habitable Earth-mass planets and moons can exist in the 

^ . Kepler-16 system, known to host a Saturn-mass planet around a stellar binary, 

by investigating their orbital stability in the standard and extended habitable 
zone (HZ). We find that Earth-mass planets in satellite-like (S-type) orbits are 
possible within the standard HZ in direct vicinity of Kepler-16b, thus constituting 
habitable exomoons. However, Earth-mass planets cannot exist in planetary- 
like (P-type) orbits around the two stellar components within the standard HZ. 
Yet, P-type Earth- mass planets can exist superior to the Saturnian planet in the 
extended HZ pertaining to considerably enhanced back-warming in the planetary 
atmosphere if facihtated. We briefiy discuss the potential detectability of such 
habitable Earth-mass moons and planets positioned in satellite and planetary 
orbits, respectively. The range of inferior and superior P-type orbits in the HZ 
is between 0.657 to 0.71 AU and 0.95 to 1.02 AU, respectively. 

Subject headings: astrobiology — binaries: general — celestial mechanics — 
planetary systems — stars: individual (Kepler-16) 



Introduction 



Kepler-16 constit utes a remarkable b inary system contai ning a circumb i nary e xtra-solar 
planet as reported by iDoyle et al.l ( 120111 ) and suggested by ISlawson et al.l (120 llh . The ex- 
oplanet was observed by the NASA Kepler spacecraft using the planetary transit method, 
which greatly enhances the confidence in the reality of the planet and provides unusually 
precise information about its physical parameters. There is a signif icant previous array 



of observations of planets in binary and multiple stellar systems (e.g., iPatience et al.l 12002 
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Eggenberger et alJbooi I2OO7I : iBonavita fc Desideral[2007l : iMugrauer k. Neuhauserll2009l : iRaghavan et al, 
2OIOI ). which also includes a small group of proposed circumbinary planets in contact binaries 
(with little chance for providing habitable environments) outside the comrnon en velopes of 
their Roche lobes Jlee et aLlboOol : loian et al.lboiol : ISeuermann et al.lboioi boilj ). 



The Kepler-16 system consists of two stars, i.e., the primary with a mass of Mi 
0.69 Mq and the secondary with M2 — 0.20 Mq, and a giant planet with properties compa- 
rable to Saturn. The orbit of the planet is almost circular and encompasses both stars. It 
takes nearly 229 days for the planet to complete one orbit. T he orbits of all thr ee objects 
are almost precisely confined to one plane (i.e., within 0.5°). iDoyle et al.l (I2OIII ) provided 
detailed information about the orbital and physical parameters of the system; however, they 
did not address issues concerning the possible habitability of the system. Obviously, it is 
the ultimate quest of the Kepler Mission to discover Earth-mass planets and moons located 
in the habitable zones (HZs) of their host stars. Therefore, it is the main objective of this 
paper to demonstrate that habitable Earth-mass planets and exomoons in stable orbits are, 
in principle, possible in the Kepler-16 system. Moreover, such objects can potentially be 
detected by the currently operating Kepler mission. This paper is structured as follows: In 
Sect. 2, we outline our theoretical approach by commenting on the stellar HZ as well as the 
employed numerical methods and considered system configurations. In Sect. 3, we describe 
our results and discussion. Our conclusions are given in Sect. 4. 



2. Theoretical Approach 

2.1. Standard and Extended Habitable Zones 

The Kepler-16 system contains two closely orbiting stars (ob = 0.22431 AU) with the 
primary (Kepler-16A) producing a substantially larger amount of photometric flux than the 
secondary (Kepler-16B) (i.e., F^/Fp^ = 0.01555). This allows us to calculate the size of the 
HZ in this system by solely taking into account the radiation of the pri mary. We compute 



the bo undaries of the HZ in this syst em by using t h e fitti ng formulas of lUnderwood et al. 



( 20031 ) based on the previous work by lKasting et al. Jl993h I also using a corrected value for 
the solar effective temperature), which are categorized by decisive atmospheric conditions of 
the Earth- mass test planet. Appropriate definitions for the inner and outer boundary of the 
stellar HZ are based on the runaway and maximum greenhouse effect, respectively, for the 
planetary atmosphere; this standard HZ is found to extend from 0.36 to .71 AU for Kepler- 
16A. This result proves to be also consistent with the more recent work by lSelsis et al.l (120071 ) 
for the case of recent Venus to early Mars-type conditions. In addition, we use the work by 
Mischna et al.l (|2000[ ) to calculate the so-called extended HZ, which requires a more extreme 
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planetary atmosphere with significantly enha nced back- warming; its outer boundary extends 
out to about 2.0 AU in the Solar System (IMischna et al.l I2OOOI ) and, correspondingly, to 
1.02 AU in Kepler-16. With the giant planet (Kepler-16b) positioned at 0.7048 ±0.0011 AU, 
it is found to be located very close to the outer edge of the standard HZ, but well within the 
extended HZ (see Fig. [1]). 



2.2. Numerical Methods and System Configurations 



The main aim of this paper is to investigate numerically the orbital stability of an 
Earth- mass object (i.e., exoplanet or exomoon) in both the standard and extended HZ of 
the Kepler-16 system. Our numerical methods are based on the Wisdom-Ho lman mapping 
technique as well as the Gragg-Burlisch-Stoer algorithm ( iGrazier et al.lll996l ). We integrate 
the resulting equations of motion forward in time for 1 million years using a fixed/initial 
(WH/GBS) time step of 10~^. The relative error in energy is calculated to determine when 
the integration methods fail and the onset of orbital instability occu rs. A secorid che ck for 
stability is performed using the method of Lyapunov exponents (see IWolf et al.lll985[ ) with 
a special emphasis on the maximum Lyapunov exponent (MLE). Our numerical simulations 
are designated as stable when they terminate with a relative energy error smaller than 10~^ 
and possess a MLE that is asymptotically approaching zero (lOuarles et al.ll201l[ ). 



By adding an Earth-mass object we consider the Kepler-16 system as a 4-body system; 
therefore, the Earth-mass object can possibly exist in 6 different orbital configurations. The 
object can move outside the orbits of the two more massive objects in a planetary (P)- 
type orbit or around any one of the massive components in a satellite (S)-type orbit. The 
putative Earth-type object could possibly exist in one of the following classes: 3 S-type 
orbits, 2 P-type orbits, or 1 Trojan exomoon. The S-type configuration would correspond 
to orbits around Kepler-16A, Kepler-16B, or Kepler-16b, i.e., the Saturnian planet. Due 
to the general definition of a moon, any S-type orbit revolving around an exoplanet would 
inherently constitute an exomoon. The P-type orbits would be classified as being either 
inferior or superior to the Saturnian planet. The Trojan exomoon orbit could exist at either 
equilateral equilibrium point, L4 or L5. 

The nu r nerica l setup of our simulations is based on the system parameters presented by 
Doyle et al.l ( 1201 ll ). see Table (H as well as adequate initial conditions. We chose an initial 
configuration with the more massive stellar component (Kepler-16A) near the center of our 
barycentric coordinate system. The less massive stellar component (Kepler-16B) is initialized 
to the left of Kepler-16A at the apastron starting position so that the initial separation of 
the stars is ab(l + eb) = 0.260 AU. Kepler-16b is initialized to the right of the primary at the 
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apastron starting position using the parameters (ap,ep) given by iDoyle et al.l (1201 ll ). The 
Earth-mass object is initiahzed to the right of the primary at the apastron starting position. 
All bodies in this system are given initial velocities in the counter-clockwise direction relative 
to the center of mass using the known eccentricities of the respective bodies. The initial 
conditions of the test planet are chosen with respect to the initial starting distance oq and 
eccentricity cq. 

In our simulations we refrain from considering S-type orbits around the stellar compo- 
nents as well as P-type orbits with initial semi-major axes less than the inner boundary of 
the standard and extended HZ. Our aim is to find stable orbits for the test Earth-mass planet 
within the standard and extended HZ. Therefore, in our computations the parameter oq is 
selected to range from 0.36 to 1.02 AU in increments of 0.001 AU, allowing us to investigate 
possible P-type planetary orbits as well as S-type orbits for the exomoon; furthermore, the 
parameter eo is selected to range from 0.0 to 0.5 in increments of 0.01. We also investi- 
gate the case of Trojan exomoons with parameters (oq, cq) equal to the initial parameters of 
Kepler- 16b, where the exomoon is placed in a position that corresponds to a point preceding 
Kepler- 16b by 60°. 

In addition to our full simulations for the 4-body system, we also give an estimate of 
the outer limit of the S-type orbital stability boundary and th e inner limit of the P-typ e 



orbital stability boundary using the statistical fitting formulas of iHolman fc Wiegertl (119991 ). 
These fitting formulas have been deduced using a range of mass ratios, distance ratios, and 
eccentricities of binaries regarding the elliptical restricted 3-body problem. It is important 
to note that the eccentricity of the test mass is neglected in this study. However, in our 
application to the elliptical restricted 4-body problem, these fitting formulas will only allow 
an estimate for test masses of low eccentricity. Since the perturbations due to the giant 
planet will be small compared to those of the stellar binary components, we implement these 
formulas neglecting the presence of this planet. This provides a good approximation of the 
conditions for stability for the S-type orbits around either of the stars, as well as the P-type 
orbits in close proximity to the inner boundary of the standard HZ. The obtained estimates 
are then used to guide our numerical study. 



Results and Discussion 



Although we refrained from considering S-type orbits around the stellar components 
(see Sect. 12. 2p . it is possible to estimate the st ability limit of S-type orb its near the stellar 



components using a statistical fitting formula (IHolman fc Wiegert 



19991 ). Our calculations 



demonstrate that an Earth-mass planet cannot exist farther than 0.0675 ± 0.0039 AU from 
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the stellar primary (Kepler-16A) for an S-type orbit due to the perturbations initiated by 
the stellar secondary (Kepler-16B). This shows that S-type orbits of a habitable Earth-mass 
planet around either of the stars must be excluded because the stability limit is well within 
the inner boundary of both the standard and extended HZ (see Sect. I2.ip . Concerning the 
P-type orbits inferior to the orbit of the giant planet, our results demonstrate that such 
orbits are unstable if the semi-major axis is smaller than 0.657 ± 0.011 AU with respect to 
the stellar primary. However, inferior P-type orbits are still possible if the test Earth-mass 
planet has a sufficient eccentricity, allowing the giant planet Kepler- 16b to capture it as an 
exomoon. 

Since all bodies in our simulations are initialized at their respective apastron starting 
position, the test planet is also given the appropriate velocity to allow for capture (see 
Fig. [2^). This leaves the principal possibility of a habitable Earth- mass planet in a P-type 
orbit between 0.657 and 0.71 AU, which means that the planet would be located within 
the standard HZ. However, short-term secular changes for such orbits allow the giant planet 
to transfer the Earth-mass planet to an orbit outside the standard HZ within 1,000 years 
(see Fig. [T^), implying that no stable P-type orbits for habitable Earth- mass planets exist 
inferior to the giant planet. The only stable P-type orbits for Earth-mass planets are those 
located superior to the giant planet. Our results show that these orbits become stable once 
their semi-major axes are 0.95 AU or higher, which places them within the extended HZ (see 
Fig. [lb). 

The last class of orbits include those that could result in an habitable exomoon in 
either an S-type or Trojan configuration. A stable S-type orbit for such an exomoon can 
be achieved through two separate scenarios. The first scenario is based on the assumption 
that the possible exomoon formed together with the giant planet, ignoring migration. The 
second scenario is based on the assumption that the putative exomoon formed initially in a 
P-type orbit and was captured by the giant planet as a result of migration. This provides 
additional justification for considering the eccentricity Cq as a free parameter because the 
initial state of the capture remains unknown. 

The first scenario which involves the exomoon forming from a secondary circumplanetary 
dis k can be estimated t hroug h the concept of Hill stability. Similar to the estimations made 



by iHolman fc Wiegertl ( 1l999l ) for the elliptical 3-body problem, there have been simulations 



to determine an approximate stability lim it for exomoons orbiting extrasolar giant planets 



( iDomingos et al.ll2006l : iKalteneggerl 120101 ) . These previous general estimations can be used 
to guide the investigation of this case. Simulations of this scenario were performed by 
incrementing the parameter cq exactly as before, but the parameter oq has been incremented 
relative to the starting position of the Saturnian planet starting from 0.0001 to 0.0240 AU 
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in increments of 0.0001 AU. The stability boundaries are expected to occur near the Roche 
hmit (inner) and the estimated Hill limit for prograde motion (outer). 

The inner boundary is calculated using the given density (see Table [T]) of the Saturnian 
planet and that of Earth, i.e., 5.515 g cm~'^. The Earth mean density has been chosen because 
it best represents what would be considered as Earth-like. Alternatively, the density could 
be varied as a function of ao, but our usage of the density is limited to the calculation of the 
Roche limit about the Saturnian planet. Planetary formation calculations for this system 
should be considered in future studies to constrain the mass and densities of possible moons 
orbiting Kepler-16b, which are however outside the scope of this paper. We also calculated 
the outer boundary assuming a circular orbit (i.e., cq = 0.0) because it gives the maximum 
stability as more ecc entric orbits have been shown to decrease the stability region (e.g.. 



Domingos et al.ll2006[ ). Thus the inner and outer boundaries are estimated as 0.0004 and 



0.0168 AU, respectively. 

Through our numerical simulation we have determined the stability boundaries to be 
commensurate with the previous estimates. The stability analyses of these orbits have un- 
dergone similar scrutiny in the relative error in energy and maximum Lyapunov exponent. 
Although we investigated the possibility of eccentric orbits in this scenario, our results show 
that the stability decreases dramatically as eccentricity is increased so that most of the stable 
orbits discovered were nearly circular (see Fig. [3^). We point out that this result is also in 
accordance with the known parameters of the most massive satellites of Jupiter and Saturn. 

The second scenario which involves the exomoon attained through capture may be 
much less probable. However, we determined that there are P-type orbits, where the proper 
eccentricity can lead to capture into S-type orbits (see Fig. |2^). The results presented in this 
figure were obtained with ao = 0.619 AU and the initial starting distance Xq = 0.699 AU, 
which is superior to the estimated stability limit of 0.657 AU. This places the test planet 
within the Hill sphere, corresponding to a Hill radius of 0.034 AU from the giant planet, 
which shows that the influence of the giant planet ranges radially from 0.67 to 0.73 AU with 
respect to the center of mass. A general trend is that the circular P-type orbits inferior to 
the giant planet are generally too fast for capture; hence, we considered eccentric orbits. 
Moreover, the P-type orbits superior to the giant planet are generally too slow for capture 
resulting in either collision with the giant planet or ejection from the system. 

If migration occurred in this system during its formation, this would be the preferred 
scenario; it wo uld entail a number of ways for a putative exomoon to avoid a chaotic orbit 



( iKippingI I2OIII ) . The cons traints on natural sa tellite formation in the Solar System were 



previously investigated by (ICanup fc Wardll2006l ). who also determined the available mass in 



the circumplanetary disk. However, the latter may not be applicable to the more exotic case 
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of the Kepler- 16 system since the Saturn- hke planet exists in close proximity to the so-called 
snow line. This makes both scenarios to appear equally likely since the migration distance 
may have been small. Further planetary formation investigations need to be performed to 
determine a more unique solution. 

Finally, we considered the possibility of a Trojan exomoon. The stability of this possible 
exomoon is ensured by the fact that the mass ratio (/i = 1112 /{nii + 1112) with nii = Mi + M2 
and 7712 = Mp) can be calculated and be used as a stability condition for the approximate 3- 
body problem. This approximation proves to be valid because the proposed Trojan exomoon 
would always reside at a semi-major axis commensurate with Kepler-16b at the appropriate 
distance leading or trailing Kepler-16b (see Sect. 2.2). Considering the case where the binary 
system can act like a point source at the center of mass, we then find the mass ratio, 
/ i = 0.000357, which is much less than the critical value for Trojans, fiQ = 0.03852 . . . 



( jSzebehelyi 119671 ) . Our results obtained for this configuration show that stable Trojans can 
exist even if the perturbations by Kepler-16B are not completely negligible (see Fig. Wp)- 
The influence of Kepler- 16B can transfer the proposed Trojan exomoon from its equilibrium 
point but it is insufficient to create an instability; instead, the Trojan exomoon would precess 
along the orbit of Kepler- 16b. 



4. Conclusions 

By performing numerical simulations we explored the possibility of a habitable Earth- 
mass object in the Kepler-16 system. Although it is beyond the scope of this paper to 
determine whether such an object exists, we are able to provide information through our 
orbital stability analyses where to search for it in the realm of further observations. We 
considered the standard and extended HZ of the system and investigated 6 different orbital 
configurations. The obtained results show that S-type planetary orbits about either of the 
stellar components as well as P-type planetary orbits inferior to Kepler-16b's orbit exhibit 
short-term orbital instabilities. These instabilities would inhibit any form of habitability due 
to the ejection, collision, or the occurrence of highly eccentric orbits for any possible object 
(if it had formed). Our main result about habitable Earth-mass planets in this system is 
that the only stable P-type orbits for such planets are those located superior to Kepler-16b's 
orbit. Specifically, these orbits become stable once the semi-major axis is 0.95 AU or higher, 
which places them within the extended HZ. The range of habitable inferior and superior 
P-type orbits in the HZ is between 0.657 to 0.71 AU and 0.95 to 1.02 AU, respectively. 

Our numerical simulations of S-type orbits for an Earth-mass moon captured by the 
giant planet or formed through a circumplanetary disk show that such orbits are stable and 
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moreover located in the standard HZ. A highly relevant aspect of the existence of an exomoon 
orbiting Kepler- 16b is that it may help explain the variation b etween the observ ed eccen- 
tricity and its previously reported value from numerical study (jPoyle et al.ll201ll ). Similar 
to the effects of Earth's Moon, the exomoon could exert possible tidal forces, through a sta- 
bilizing torque and angular momentum exchanges, freezing the eccentricity of Kepler-16b at 
the observationally determined value of 0.0069. As the observations may just be a snapshot 
in time, other hypotheses concerning the Kozai Mechanism, driving of eccentricity through 
planet-planet interactions, or possible exomoons, can be considered equally likely and can- 
not be uniquely id e ntified until further invest igations of the system have been performed 
( Laskar et al. 1993 : Barnes &: Greenbergl 2008). Since the exomoon would lie on the outer 
edge of the HZ, tida l heating could also pl ay a role in increasing the prospects of habitability 
( iBarnes et al.l l2009l : [Dressing et al.ll2010l ). We also considered the possibility of a Trojan 
exomoon present at one of the equilateral Lagrange points. This is an intriguing case as 
the exomoon would precess and display a variety of different orbits about L4 or L5, which 
are all located within the standard HZ. Hence, the existence of habitable exomoons around 
Kepler-16b is an exciting scenario for facilitating habitability in the Kepler-16 system. 

Finally, we want to point out that an Earth-mass planet or moon can pot entially be de 



tected by the current Kepler mission. Our suggestion is based on the work by lKipping et al. 



( 120091 ) who provided evidence that Kepler's instruments should indeed be capable of detect- 
ing such objects, specifically, the possible exomoon. In general, Kepler should be able to 
de tect exomoons a s sma ll as 0.2 if observers look for the transit timing effects discussed 
by iKipping et al.l ( 120091 ). In addition, the realm of circumbinary planets similar to Kepler- 
16b may prove to be an adequate region for discovering exomoons as it would be possible to 
obtain the necessary statistical constraints of detection. The 220-day period the Saturnian 
planet makes it possible to observe 6 transits within 3.6 years of observation, which is well 
within the extended mission life time of the Kepler spacecraft. It is also noteworthy that 
the distance to the Kepler-16 is only about 61 parsec, which places the system well within 
the distance range (0.68 - 386 pc) for detecting Earth-mass planets and exomoons around 
K and M-type stars. 



This work has been supported by the U.S. Department of Education under GAANN 
Grant No. P200A090284 (B. Q.), the Alexander von Humboldt Foundation (Z. E. M.), and 
the SETI institute (M. C.). Moreover, it has been supported in part by NASA through the 
American Astronomical Society's Small Research Grant Program (M. C.). We also would 
like to thank R. Heller (the referee) for valuable comments. 



- 9 - 



REFERENCES 

Barnes, R., & Greenberg, R. 2008, in lAU Symp. 249, Exoplanets: Detection, Formation 
and Dynamics, ed. Y.-S. Sun, S. Ferraz-Mello, & J.-L. Zhou (Cambridge: Cambridge 
University Press), 469 

Barnes, R., Jackson, B., Greenberg, R., & Raymond, S. N. 2009, ApJ, 700, L30 

Beuermann, K., et al. 2010, A&A, 521, L60 

Beuermann, K., et al. 2011, A&A, 526, A53 

Bonavita, M., & Desidera, S. 2007, A&A, 468, 721 

Canup, R. M., & Ward, W. R. 2006, Nature, 441, 834 

Domingos, R. C, Winter, O. C, & Yokoyama, T. 2006, MNRAS, 373, 1227 
Doyle, L. R., et al. 2011, Science, 333, 1602 

Dressing, C. D., Spiegel, D. S., Scharf, C. A., Menou, K., & Raymond, S. N. 2010, ApJ, 721, 
1295 

Eggenberger, A., Udry, S., & Mayor, M. 2004, A&A, 417, 353 

Eggenberger, A., Udry, S., Chauvin, G., Beuzit, J.-L., Lagrange, A.-M., Segransan, D., & 
Mayor, M. 2007, A&A, 474, 273 

Grazier, K. R., Newman, W. I., Varadi, F., Goldstein, D. J., & Kaula, W. M. 1996, DDA 
Meeting No. 27, BAAS, 28, 1181 

Holman, M. J., & Wiegert, R A. 1999, AJ, 117, 621 

Kaltenegger, L. 2010, ApJ, 712, L125 

Kasting, J. F., Whitmire, D. R, & Reynolds, R. T. 1993, Icarus, 101, 108 
Kipping, D. M. 2011, Ph.D. thesis. University College London (larXiv:1105.3189|l 
Kipping, D. M., Fossey, S. J., & Campanella, G. 2009, MNRAS, 400, 398 
Laskar, J., Joutel, F., & Robutel, P. 1993, Nature, 361, 615 

Lee, J. W., Kim, S.-L., Kim, C.-H., Koch, R. H., Lee, C.-U., Kim, H.-L, & Park, J.-H. 2009, 
AJ, 137, 3181 



Mischna, M. A., Kasting, J. F., Pavlov, A., & Preedman, R. 2000, Icarus, 145, 546 
Mugrauer, M., & Neuhauser, R. 2009, A&A, 494, 373 
Patience, J., et al. 2002, ApJ, 581, 654 
Raghavan, D., et al. 2010, ApJS, 190, 1 

Selsis, P., Kasting, J. P., Levrard, B., Paillet, J., Ribas, I., & Delfosse, X. 2007, A&A, 476, 
1373 

Slawson, R. W., et al. 2011, AJ, 142, 160 

Szebehely, V. 1967, Theory of Orbits (New York and London: Academic Press) 
Qian, S.-B., Liao, W.-R, Zhu, L.-Y., & Dai, Z.-B. 2010, ApJ, 2010, L708 
Quarles, B., Eberlc, J., Musielak, Z. E., & Cuntz, M. 2011, A&A, 533, A2 
Underwood, D. R., Jones, B. W., & Sleep, P. N. 2003, Int. J. AstrobioL, 2, 289 
Wolf, A., Swift, J. B., Swinney, H. L., & Vastano, J. A. 1985, Physica D, 16, 285 



This preprint was prepared with the AAS macros v5.0. 



- 11 - 



Table 1: Stellar and Planetary Parameters of Kepler- 16 



Parameter 


Value" 


Distance (pc) 


~ 61 


Fb/Fa 


0.01555 ± 0.0001 


Ml (Mq) 


0.6897 ± 0.0035 


M2 (Mq) 


0.20255 ± 0.00066 


Teff,l (K) 


4450 ± 150 


Ri (Rq) 


0.6489 ± 0.003 


^b(d) 


41.079220 ± 0.000078 


ab (AU) 


0.22431 ± 0.00035 


eb 


0.15944 ± 0.00061 


Mp (Mj) 


0.333 ± 0.016 


ap (AU) 


0.7048 ± 0.0011 




0.0069 ± 0.001 


Pp (g cm-3) 


0.964 ± 0.047 



Note. — "Data as provided by 



Doyle et al 



(|201l[ ). All parameters have their usual meaning. 
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Fig. 1. — (a) Depiction of an unstable P-type Earth-mass planet (black) with an initial 
semi-major axis of = 0.504 AU and an initial eccentricity of Cq = 0.06 to give a starting 
position at apastron of xq = 0.534 AU; the stars are given as A (red) and B (blue), (b) 
Depiction of a stable P-type Earth- mass planet (black) with an initial semi- major axis of 
ao = 0.951 AU and an initial eccentricity of eo = 0.03 to give a starting position at apastron 
of Xq = 0.980 AU. The dark gray region represents the standard habitable zone (HZ) and 
the light gray region represents the extended habitable zone (EHZ). The agreement between 
the orbit of the giant planet Kepler-16b (purple) and the outer edge of the standard HZ is 
coincidental. 
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Fig. 2. — (a) Depiction of an S-type captured Earth-mass exomoon (black) with an initial 
semi-major axis of qq = 0.619 AU and an initial eccentricity of eo = 0.13 to give a starting 
position at apastron of xq = 0.699 AU; the stars are given as A (red) and B (blue), (b) 
Depiction of a possible Trojan exomoon in a rotating reference frame (black) with an initial 
semi-major axis of = 0.7048 AU and an initial eccentricity of Cq = 0.0069 to give a 
starting position at apastron of xq = 0.710 AU. The dark gray region represents the standard 
habitable zone (HZ) and the light gray region represents the extended habitable zone (EHZ). 
The agreement between the orbit of the giant planet Kepler- 16b (purple) and the outer edge 
of the standard HZ is coincidental. 
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Fig. 3. — (a) Depiction of a stable S-type coformed Earth-mass exomoon (black) with an 
initial semi-major axis of qq = 0.715 AU and an initial eccentricity of cq = 0.0; the stars 
are given as A (red) and B (blue), (b) Depiction of an unstable S-type coformed Earth-mass 
exomoon (black) with an initial semi-major axis of Oq = 0.721 AU and an initial eccentricity 
of eo = 0.0. The dark gray region represents the standard habitable zone (HZ) and the light 
gray region represents the extended habitable zone (EHZ). The agreement between the orbit 
of the giant planet Kepler- 16b (purple) and the outer edge of the standard HZ is coincidental. 



